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Abstract 

This paper is devoted to studies of the accuracy of approximate equations for 
the reflection function of semi-infinite turbid media. It is assumed that the 
probability of photon absorption 6 by a particle is low. We derive quadratic 
and cubic terms of the expansion of the reflection function with respect to 
./B. Our results show that the cubic term makes a great improvement in the 
accuracy in the region of values of the similarity parameter that are not small. 


1. Introduction 


The reflection function of a semi-infinite turbid medium plays 
a central role in the radiative transfer theory [1,2]. It obeys 
the nonlinear integral equation [1]. Nowadays there is no 
problem to calculate this function numerically and codes are 
freely available via Internet [3]. 

This in no way diminishes the importance of analytical 
results. In particular, such results can be obtained in the 
case of small probabilities 6 of photon absorption by particles 
[4-9]. For instance, it follows that for the reflection function 
Ro(&, n, p) of a semi-infinite plane parallel turbid medium as 


B— 0, 


Roo(E,0, p) = RVE, n, o) — yKo(E)Ko(m), Œ) 


where R? E, n, 9) = Ræ(é, n, p) at B = 0, £ is the cosine of 
the incidence angle, 7 is the cosine of the observation angle, g 


is the azimuth angle, 
3(1 — g) 


is the similarity parameter and 


=> 


1 JE 
Í p(0) sin 0 cos 6 dé (3) 
2 Jo 


is the asymmetry parameter, which describes the degree of 
elongation of the phase function p(@) in the forward direction. 
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The function Ko(&) describes the angular distribution 
of light escaping from a turbid semi-infinite plane-parallel 
nonabsorbing layer with sources of radiation, located at infinity 
[6]. This is a smooth function of the escape angle due to the 
randomization of photon directions during their propagation 
from infinity to the border of a semi-infinite nonabsorbing 
medium. Actually Ko(&) practically does not depend on the 
angular distribution of light in a single scattering event p(0) 
(at least for € > 0.2). It follows, approximately, that [5] 


Ko(&) = 3(1 +28). (4) 


We see that equation (1) allows us to relate the reflection 
function of an absorbing turbid medium as 6 — 0 with that 
for a nonabsorbing turbid layer. 

This result is of a considerable importance due to the 
fact that the function R? (€, n, 9) depends only on the phase 
function of a scattering medium and parameters €,7,@. It 
does not depend on £ by definition. Moreover, for selected 
classes of turbid media and viewing geometries, the function 
R? (E, n, o) practically does not depend on the specific form 
of the phase function as well [10, 11]. This allows us to make 
the substitution of this function by a simple universal function, 
which depends only on the geometry [11]. 

Note that equation (1) can be applied to the determination 
of the spectral dependence y(A), where à is the wavelength. 
This is of importance for the spectroscopy of turbid media [10]. 

The task of this paper is twofold. First of all, we derive 
next two terms of the expansion, given by equation (1), which 
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are proportional to y? and y? correspondingly. For this we 
will use well-known expansions for a plane albedo of a light 
scattering medium [7,8]. Secondly, equation (1) and new 
expansions are compared with radiative transfer calculations 
using the code described by Mishchenko et al . 


2. Theory 


Clearly, in the general case for arbitrary values of y, we have, 
instead of equation (1), 


Rol, n, g) = any”, (5) 
n=0 


where, comparing equations (1) and (5), we have: aọ = 
R? (E, n, g) anda, = —Ko(&)Ko(n). Coefficients a, atn > 1 
are unknown. 

Our task is to find values of a2 and a3. For this, we will 
use the well-known expansion for the integral 


1 20 1 
r) = zf a f ndn RælĖ, n, p), (6) 
T Jo 0 


whichis called the plane albedo. It gives us the total reflectance 
of a turbid layer at a given light incidence angle Vo = arccos £. 
Namely, it follows at small y that [7, 8] 


r€) = Da". (7) 
where 
by = 1, bı = —Ko(é), 
by = 3(pKo(€) + 54 Vo(é)), (8) 
bs = —ġ (u — vWo(E)) Kol) 
and 
u=1—3g—(e-#)q, v = (e — hq — åp’, 
p= af KEE de, q= E, 
0 l-h 
e= of Ko(E)E? dé, (9) 


h an (0) sin@ cos? 6 dé : 
=- sin@ cos — —, 
Ady? 2 


1 Qn 1 
we=e-— f ay | R? Œ, n, on? dn, 
T Jo 0 


Ko) 
2Ko(E) fy Kamn dn 


Here K(&) is the coefficient in the expansion of the escape 
function, given by 


Wo(§) = 


KE) =) K, E)". 


n=0 


(10) 
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It follows that [8] 


Kı) = -p Ko), (11) 


where Ko(&) at € > 0.2 is given by equation (4). 

The angular dependence K>(&), unlike the case of K; (£), 
differs from Ko(&). One can find functions Vo(é) and Wo(&), 
which define K2(&) (see equations (9)) in figures 1 and 2. 
Symbols represent the data of Minin [8] for cloudy media and 
solid lines correspond to our approximate equations, obtained 
by fitting Minin’s [8] numerical data. They are 


Vo(E) = A187 — 36 — 2], 


i ip 20542 (12) 
WoE) = 3+ 3E + RE 


Importantly, equations (12) can be applied to media other 
than clouds. This is possible due to their low sensitivity to 
the detailed microstructure of a turbid medium in question 
[8]. Integrals for p and ¢ can be found analytically, using 
equation (4), which is also not sensitive to the microstructure 
of a disperse medium. The simple integration gives us 


10 -u 
E 


(13) 


0.47 


0.2 F 


0.0 F 


Figure 1. Function Vo(&) calculated with equation (12) as compared 
to the numerical data of Minin [8] for cloudy media. 


Figure 2. Function Wo(&), calculated with equation (12) as 
compared to the numerical data of [8] for cloudy media. 


Accuracy of approximations for the reflection function 


Thus, we have (see equations (9)) 


u=1+2q—- 3g, (14) 
We see that parameters u and v in equation (8) depend solely 
on values of g, h and the ratio q = (1 — g)/(1 — A), which 
are determined by the phase function of a scattering medium 
in question. 

In particular, it follows that for characteristic sizes 
of water droplets in terrestrial clouds: g = 0.84—0.87, 
h = 0.76-0.80, q = 0.65-0.67. These results were obtained 
using Mie calculations for the Y water droplet size distribution 
Yaf exp(—9a /aet), where Y is the normalization constant and 
er is the effective radius, defined as the ratio of the third to the 
second moments of the droplet size distribution. The droplet 
effective radius aef was changed in the range 4-20 um, which 
is characteristic for water clouds. The wavelength à was equal 
to 0.65um. 

One can see that parameters g, h and q are not particularly 
sensitive to the microstructure of a cloudy medium. So for 
rapid estimations, one can use fixed values of g, h and q. For 
instance, using values of g = 0.84, h = 0.76 and g = 2/3, 
one obtains, for b) and b; in equation (7), 


by = aKo(§) + y Vol), b3 = —(@+ y Wo(E)) Ko). 
(15) 


Here we have (see equations (8), (13) and (14)) 


for water clouds. The values of a, y, a and y for other types 
of disperse media can be easily found if integrals g and h (see 
equations (3) and (9)) are known. 

So instead of equation (7), it follows that 


r(&) = 1 — yKo(&)+y*(aKo(&) + y Vo(&)) 
—y'(&@ + PWE) KoE). 


Let us generalize equation (16) and consider the case 
of the reflection function R.o(€,7, ) in the framework of 
the same approach. Our derivation is not mathematically 
rigorous, but rather physically based. Its accuracy, however, 
appeared to be high as compared to the numerical solution 
of the radiative transfer equation (see next section). For this 
we will use equations (6) and (16), the reciprocity principle 
(Ro(—,7,¢9) = Ron, &,y)) and following exact integral 
relations [8]: 


(16) 


1 
2 f KoE)WE)E dé = 1, a7) 
1 
af KoE dé = 1, (18) 
0 
1 
[ Vo(é)é dé = 0, (19) 
1 Qn 1 
zf a f ndnR9,(,7, 9) =1. (20) 
T Jo 0 


Then instead of equation (16), we have 


Roo(&. n, ) = R(E, n, P) — y KoE) Kon) 
+y” le Ko(§)Ko(n) + y(Vo(é) + Vo(n))] 
—y*L(@ + ¥ Wo(E)Wo(n)) Ko(E)Ko(n)I. 

Equation (16) follows from equation (21) using the 

definition (6) and integral relations (17)-(20). Due to the 


reciprocity principle, we can have in equation (21) only 
combinations 


(21) 


Ko(&) + Ko(n), Kol) Ko(n), 
Vo(&) + Vo(n), Vo(&)Vo(n), 
Wo(&) + Wo(n), Wo(&) Wo(n). 


Combinations Ko(&) + Ko(7), Vo(&) Vo(n) and Wo(&) + Wo(7) 
do not enter equation (21) because they are not consistent with 
equations (6) and (16). Note that equation (21) allows us to 
find the spherical albedo 


2 T 1 1 
T Jo 0 0 
Namely, it follows that 


r, = l—y+tcy* — dy’, (23) 

where c = a ~ 0.54, d = & +y ~ 0.15 for clouds or without 

substitution for values of p, q, u, v for arbitrary disperse media: 

c= Èp, d= (u — v). (24) 

Equations (23) and (24) were earlier derived by Yanovitskij 

[7]. This confirms our calculations. Note that equation (23) 

reminds us of the following well-known expansion: 

e? = 1— y+A1y — Ay? +s, (25) 

where A, = 1/2, Ap = 1/6. Indeed we have, approximately, 
c © A\,d X A and most importantly, 


=e? (26) 


This approximate phenomenological equation is known for 
a long time [4]. One can see that this simple formula 
approximately accounts for up to the third term in the 


expansion 
ee] 
To = a Nay”. 
n=0 

As we will see, it can be used even beyond the range of vadility 
of the exact asymptotic result (23) (particularly, at y ~ 1). Note 
that equation (26) can be generalized to describe the reflection 
function Ro(é, n, o) [4,9]: 


Roo (E,n, p) = RQ (E, n, p) exp(—u(E,n,g)y), 27) 
where KoE) Kolm) 
0 on 

Sh; = —_.. 28 

AE = RE E) =a 
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3. Accuracy of approximations 


Let us study now the accuracy of equations (1), (21), (23), 
(26) and (27). For this we will use the numerical solution 
of the Ambarzumian’s nonlinear integral equation for the 
reflection function of a semi-infinite random medium [3]. This 
equation is valid for arbitrary values of the probability of 
photon absorption. The results of exact calculations together 
with data obtained from equations (23) and (26) are presented 
in figure 3 as functions of the similarity parameter y. One can 
find errors A of these equations in figure 4. 

Approximate equations and maximal values of y and 
@ = 1 — B, which allow to use them with accuracy better 
than 5%, are presented in table 1. We find also that if the 
tolerable error is 10%, then maximal values of y are equal to 
0.45, 0.65, 1.0 and 1.4, depending on the approximation used 
(see table 1). 

Note that taking into account the cubic term allows us to 
increase the value of y™* by 2.5 times (see table 1) as compared 
to the linear approximation (1). This is of importance for 
many natural media, including terrestrial and venerian clouds. 
The best accuracy is, however, due to the simple exponential 
solution (26) [4]. 

The reflection function R(1,0.6,0) is presented in 
figure 5. Again we see that the exponential solution (27) has 
larger range of applicability than equation (21). It is simpler 
in structure also. 

It should be stressed that data in figures 3 and 5 were 
obtained for different values of aef and, correspondingly, for 
different phase functions. They perfectly lie on a smooth 


Egeru bao eee ale eer 


0. 
0.0 05 10 15 2.0 25 3.0 3.5 40 45 5.0 
y 


Figure 3. Dependence r, (y) according to equation (26) (line 1), 
(23) (line 2—the account only for the linear term, line 3—the 
account for the quadratic and linear terms, line 4—the account for 
the linear, quadratic and cubic terms) and exact radiative transfer 
calculations (5) for cloudy media at A = 1.549 um, aef = 4, 6 and 
16 um, and the refractive index of droplets n = 1.3109 — ix, where 
x = 1075, 1074 and 107°. 
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Figure 4. Errors, correspondent to data, in figure 3 (1—linear term, 
2—linear and quadratic terms, 3—linear, quadratic and cubic terms, 
4—exponent (26)). 


Table 1. The maximal values of y and wo, which make the accuracy 
of equations better than 5%. The value of g was equal to 0.85. 


Roo ymax ay 

l-y 0.3 0.9975 
1 — y +cy? 0.5 0.9930 
1—y+cy? — dy? 0.75 0.9842 
ey 1 0.9719 
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Figure 5. Dependence of the reflection function Ra (1, 0.6, 0) on 
the similarity parameter y, obtained from the exact radiative transfer 
calculations (1), equation (21) (2—the account only for the linear 
term, 3—the account for the linear and quadratic terms, 4—the 
account for the linear, quadratic and cubic terms) and equation (27) 
(stars) for cloudy media at A = 1.549 wm, aef = 4, 6 and 16 um, 
and the refractive index of droplets n = 1.3109 — ix, where 

x = 1075, 1074 and 107°. 
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Figure 6. (a) Error of approximations, given by equation (21) 
(1—linear term, 2—linear and quadratic terms, 3—linear, quadratic 
and cubic terms) and equation (27) (triangles). The input parameters 
are the same as in figure 5, but aef = 6 um, x = 1074. (b) Same as 
in figure 6(a), but dep = 4 um, x = 5 x 1074 (1—linear and 
quadratic terms, 2—linear, quadratic and cubic terms, 3— 

equation (27)). (c) Same as in figure 6(a), but aef = 6 um, 

x = 5 x 1074 (1—linear, quadratic and cubic terms, 

2—equation (27)). 


curve as plotted against y. This suggests that media with 
different values of wo and g but the same y have similar 
radiative characteristics. This is the numerical confirmation 
of the similarity principle [6]. 


Table 2. Values of wo, g and y for different aef and x as obtained 
from the Mie theory at A = 1.549 um, n = 1.3109. 


Aef X w0 g y 

6 1074 0.9953 0.8210 0.37 
4 5 x 1074 0.9850 0.7912 0.51 
6 5 x 1074 0.9774 0.8253 0.83 


The errors of approximations (21) and (27) for the 
reflection function at nadir observation and different angles 
of incidence are presented in figures 6(a)-(c). They were 
obtained at A = 1.549 um, ag = 6m, n = 1.3109, 
x = 10-4 (figure 6(a)), dep = 4 um, x = 5x 1074 (figure 6(b)) 
and dep = 6 um, x = 5 x 1074 (figure 6(c)). The values wọ, g 
and y for all three cases were calculated with the Mie theory. 
They are presented in table 2. 

We see that the error of the approximation depends on 
the incidence angle. It drastically increases for grazing light 
incidence. Note that equation (21) provides the best accuracy 
at small y. On the other hand, equation (27) has wider 
applicability range (see figure 6(c)) and can be applied even 
for relatively large values of y, where the total reflectance of 
light from a layer is low. 


4. Conclusion 


We have studied here the accuracy of different equations 
for the fundamental function of the radiative transfer theory, 
Ro (&, 7, o). It was found that equation (21) is more accurate 
as B — 0. However, the use of equations (26) and (27) 
is preferable if one needs to make rapid estimations or the 
interpretation of experimental data especially at y ~ 1. 
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